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Abstract. We study the attractive Bosc-Hubbard model with a tunable, on-site 

three-body constraint. It is shown that the critical behavior of the system undergoing 

a phase transition from pair-superfiuid to superfluid at integer filling depends on the 

value of the three-body repulsion. In particular, a critical exponent and the central 

charge governing the quantum phase transitions are shown to have repulsion dependent 

features. In consequence, the model extends the list of known systems violating the 

universality hypothesis. 
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1. Introduction 

Amazing experimental progress on ultra-cold atoms confined in optical lattices has 
generated an almost ideal arena for simulating and testing various Hubbard-like models 
[H [21 |3j. Since the first experimental confirmation of the quantum phase transition 
(QPT) from the Mott insulator phase (MI) to the superfluid phase (SF) in the Bose- 
Hubbard model [4], many different extensions of the model were proposed theoretically 
to study [21 IB]- Extended Hubbard models (for fermions and bosons) take into account 
processes induced by internal structure of interacting particles, long-range forces, and 
higher bands of optical lattices. All of these lead directly to additional interaction terms 
in the Hamiltonian which support a rich plethora of physical properties [5] awaiting 
experimental verification. 

Availability of such systems in lower than three dimensions promises avenues to 
investigate the nature of phase transitions in low dimensionality. Particularly, for ID 
quantum systems it opens a new route to look for and test conformally invariant systems 
with different central charges c. It is known that for systems with c < 1 the critical 
exponents are fully determined by the charge and they do not depend on the details of 
the microscopic model [6]. This is a specific manifestation of the so-called universality 
hypothesis (UH) close to quantum phase transitions (QPT). However there are models 
violating the hypothesis in the sense that their critical exponents continuously vary with 
respect to the microscopic parameters of the model. Prominent examples of such non- 
universal behavior are found in Ashkin- Teller model [3 El E] which is a specific instance 
of a compactified free bosonic model generically described by the theory Tomonaga- 
Luttinger liquid theory [10], eight- vertex model solved exactly by Baxter |llj . the 
arbitrary spin Heisenberg models [l2j , etc. From the conformal invariant point of view 
they are related to field theoretical models with central charge c > 1. 

A recent unconventional idea leading to a system described by a CFT with fractional 
central charge larger than one centered around the study of attractive ultra-cold bosons 
described by a Bose-Hubbard model with an on-site three-body constraint [TBI IT4] . 
The three-body constraint dis-allows the occupation of a site by three or more bosons 
and is equivalent to an infinite on-site three-body repulsion term. Importantly, the 
constraint inhibits the collapse of bosons in the system on to a single site and the 
resulting competition between tunneling and attractive two-body interactions leads to 
interesting quantum critical behaviour. At low tunnelings, a pair-superfluid phase (PSF) 
is supported while at high tunnelings the more conventional single particle superfluid 
phase (SF) is favoured. These two phases are separated, at unit filling, by an exotic 
Ising type quantum phase transition, albeit characterized by a U(l) x Z2 conformal 
field theory with fractional central charge c = 1 + 1/2. This model is experimentally 
feasible as the three body constraint can be enforced by rapid three-body recombination 
processes corresponding to the decay into the continuum of unbound states [13] and 
thus loss of particles from optical lattices. This has driven interest in deep studies of 
the model [151 EEl [T71 HH [191 [20] . In contrast to other systems explored with the PSF 
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phase [2H [221 [231 EH US\, here the PSF emerges through local interactions only. It was 
further shown that such a model can be mimicked by using a filling-one Mott insulator 
of spin- 1 atoms [26J. 

2. The model 

In this article we consider the attractive Bose-Hubbard model with soft three-body 
on-site repulsion. We assume that the softening of the three-body repulsion is still 
sufficiently large to prevent system from collapsing. Whilst the three-body losses 
dominate in the system, they are not instantaneous. We restrict to the one-dimensional 
lattice wherein the system of L sites (with open boundary conditions) is described by 
the extended Bose-Hubbard Hamiltonian: 

L-l jj L 

H = -J^2 {aw+i + oJ+iOj) + — Y^ ni(hi - 1) 
i=i l t=i 

W J^ 
+ — J2ni{n i -l){n i -2). (1) 

6 £i 

The hi is the annihilation operator of single boson on site i and hi = a\a,i is the local 
particle number operator. The two-body interaction U < and W is the repulsive three- 
body interaction parameter. This three-body term can be understood as an effective 
way of taking into account higher orbital physics and three-body losses. In the limit 
W — > oo the last term in the Hamiltonian is equivalent to the three-body on site hard- 
core constrain Vi(aJ) 3 = introduced and studied in detail earlier [T5l fT6| ITT ] ITS ] IT9T 120] . 
Here, we examine the properties of the tunable QPTs the system undergoes when the 
a three-body term is large but finite. While populating a lattice site with three bosons 
is not completely forbidden, the associated energy cost makes this process unfavorable. 
Thus states with sites occupied by more than two bosons affect the system properties 
mainly via virtual processes. It is worth noting, that such a model was also studied 
in the opposite case, i.e. for repulsive two-body interactions and attractive three-body 
ones, where contraction of insulting lobes was predicted [23, EBJ [29] . 

We study the QPT numerically using two approaches: Exact Diagonlization (ED) 
of small clusters and the Density Matrix Renormalization Group method (DMRG) for 
larger sizes. Using these methods, we study the variation with the repulsion parameter 
W of quantum critical properties at the mentioned QPTs. 

It is known that properties of Hubbard systems crucially depend on the relation 
between total number of particles in the system iV and the number of lattice sites L. 
Therefore, it is convenient to introduce filling factor p = N/L. One can also rewrite 
the Hamiltonian to the grand canonical ensemble by introducing chemical potential /i 
and by adding the additional term — \x J2i "hi to the Hamiltonian. In this picture the 
average number of particles in the system is determined by the chemical potential. In 
Fig. [1] the phase diagram of the hard-core constrained system (W — ¥ oo) is presented. 
It is worth emphasizing that although the Hamiltonian is mathematically correct for all 
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Figure 1. (color on-line) The phase diagram of hard-core constrained model (W — > 
oo) obtained with exact diagonalization method in full many-body basis on 8 lattice 
sites, (a) average filling p, (b) single-particle superfluid hopping cj>, (c) pair-superfluid 
hopping <F For p = and p = 2 the system always is in the insulating phase. For 
intermediate filling, depending on the tunneling, the system is in SF or PSF phase. 



fillings p < 2, it describes the physical system correctly only for low fillings. At higher 
fillings the influence of higher orbitals should be taken into account for accuracy. In the 
following, we specialize to the unit filling case p — 1, i.e. when total number of particles 
N is equal to the total number of lattice sites L. 

The hard-core constrained model with W — > oo in (p]) undergoes the phase 
transition from PSF to SF for any p < 2. In the grand canonical ensemble the SF phase 
is characterized by non vanishing order parameters (a^) ^ and (a^) ^ 0. In contrast, in 
the PSF phase the condensate fraction (<2j) vanishes. Note that the difference between 
PSF phase and insulating phases (which we do not conisder further) lies in the dimer 
order parameter (af) which in the PSF remains non vanishing, i.e. nonlocal correlations 
still exist. In the studied situation, when total number of particles is conserved, 
these two order parameters can not be found directly, since they are calculated as 
expectation values of non- conserving total number of particles operators. Therefore, 
following [25], we introduce two other correlation functions (j) = J2( a i+i a i)/(L — 1) and 
$ = X)(aI+iaI+i a tOi)/(L — 1). These two correlation functions, called hopping fields, 
describe hopping of particles and correlated pairs between neighboring sites. They might 
be viewed as a number conserving analogs of the ordinary and dimer order parameters 
respectively. 



Violation of the universality hypothesis in ultra-cold atomic systems 



3. Results 
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Figure 2. Hopping fields <fi and $ obtained with the exact diagonalization of the 
Hamiltonian ([1]) for system with N = L = 10 in two extreme cases of three-body 
repulsions. PSF hopping (thin line) is always present in the system and it dominates 
for small tunneling rates. SF hopping is negligible for small tunnelings and it domanetes 
in the system for large ones. For some intermediate tunneling amplitudes the system 
undergoes a phase transition between these two superfluid phases. 



To examine the critical behavior of the system for p = 1 for different values of three-body 
repulsion we first perform ED of the Hamiltonian for given lattice size L. In this way, for 
given tunneling J (measured in the units U = 1) we find the exact ground state of the 
system |G( J)) and its energy E G (J). The simplest way of finding properties of studied 
system is to calculate hopping fields (see Fig. [2]). It is a straightforward observation 
that depending on tunneling system exhibit nonlocal correlations. For small tunneling 
SF hopping is negligible, and PSF hopping $ dominates in the system. For larger 
tunneling amplitudes the SF hopping increases and plays a crucial role in determining 
properties of the system. Apart from some details this general picture holds also even 
for very soft three-body repulsions (Fig. |2b). Hence the QPT from PSF to SF exists 
for any reasonable value of W. 

We now use the fidelity susceptibility as a quantitative indicator of the PSF-SF 
QPT, defined as [321 EH 1331 ED] : 



x(J) 



d 2 F{J, 5) 



d5 2 



6=0 



-21im 

(5->0 



In J 7 (J, 6) 

6 2 



(2) 



where J 7 (J, 5) = \(G(J)\G(J + S))\ is the fidelity of neighbouring ground state of the 
system differing by a small change of the tunneling rate J. In the thermodynamic limit, 
i.e. when L — > oo, the singular behavior of FS is an efficient indicator of the phase 
transition of finite rank. For finite system size, the FS is a smooth function of J but 
some non monotonic behavior is still present. In the left panel of Fig. [3] the FS for two 
different three-body repulsions W and for different system sizes L are presented. One 
can find that maximal amplitude of the FS, as well as its position crucially depend on the 
system size L. This behavior of the FS suggests that for infinite size system at critical 
point fidelity J 7 will collapse to zero and its susceptibility x will diverge. Therefore, 
studied quantum phase transition is of the Landau's type and it can be understood 
with phenomenological Landau symmetry-breaking theory [321 EH EH ESI E3 E3]. This 
picture is consistent with previous results for hard-core model, since it is known that in 
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such a system PSF-SF transition is of the Ising-type at unit filling p = 1 [TBI El [151 US]. 
From this behavior of FS, by using the finite scaling method (see for example 
[301 [38]). one can extract the position of the critical tunneling J c and critical exponents 
describing divergences of system properties near the critical point. To do this 
systematically we introduce the reduced tunneling t = (J — J c )/J c as the measure 
of the distance from the critical point. Then we assume that in the neighborhood of the 
critical point the correlation length £ and the fidelity susceptibility \ diverge according 
to some power laws £ ~ \t\~ v and x ~ |*|~ 7 - Combining these two relations we find the 
connection between FS and correlation length for the infinite system size, i.e. x ~ £^l v . 
For finite size systems the situation is different since correlation length £ can not be 
larger than system size L. More precisely it means that the relation is valid only as long 
as £ <C L. When £ become comparable with the system size L then the FS approaches 
some finite value and does not change. This observation can be expressed by introducing 
an additional function Xo{ x ) which has some asymptotic value K in the limit x — > oo 
and it diverges as x 1 ^ when x — >• 0. Then one can write x = £ 7 ^Xo (£/£)• To make this 
relation independent of the correlation length, we finally introduce the scaling function 
x(x) = x~' y Xo( xV )- In this way we can write susceptibility as a function of reduced 
tunneling t and system size L 

X (t,L)=L^x(L^t). (3) 

The scaling function x is analytical and in principle is a universal function independent 
of the system size. In practice, it contains some dependence on system size far from its 
extrema. 

Equation ([3]) can be used to determine critical tunneling J c and critical exponents 
v and 7. For example, there is a very simple way of determining the relation between 
two critical exponents 7 and v. By plotting the maximal amplitude of the fidelity 
susceptibility Xmax as a function of system size L on a log-log scale one expects linear 
relation with a slope equal to j/u. Plots in the Fig. H] show that even for small 
systems the linear regression is almost perfectly satisfied for the hard constraint model 
(W — > 00) as well as for very soft three-body repulsion (W = 4|[/|). Absolute values of 
the critical exponents as well as the value of the critical tunneling J c can be found via 
the data collapse method. It follows via the following simple observation. If all critical 
parameters were known and indeed the relation ([3]) would hold then by plotting L^^x 
against L l ^t for different system sizes all curves should collapse to the universal function 
X- Moreover, this collapse occurs only for unique set of critical parameters. Therefore, 
in principle in this way all critical parameters can be determined. The example data 
collapsing for studied system in two cases W = 4\U\ and W = 00 are presented in 
right panel in Fig. [3j In Table [TJ the critical parameters for intermediate three-body 
interactions W are presented. These results show that, as expected, the position of 
the transition point from pair to normal superfluid phase depends on the three-body 
repulsion. Interestingly however, in addition we find that the critical exponent 7 varies 
with W and is thus dependent on the microscopic details of the model. This means 
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Figure 3. (color on-line) Fidelity susceptibility x as a function of tunneling amplitude 
J for two extreme values of the three-body repulsions (W — 4|J7| and W — oo) obtained 
for different system sizes L — 6,8, 10, 12. (left panel) FS is non monotonic function of 
J and its peaks are indicators of changes in the ground state of the system (quantum 
phase transitions), (right panel) the same plot with rescaled variables. Rescaling 
parameters (critical tunneling J c and critical exponents v and 7) for which all curves 
collapse to one universal curve are given on the right. The critical exponent v is 
insensitive to the strength of a three-body repulsion. In contrast, critical exponent 7 
as well as critical tunneling J c depends on W. 
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Figure 4. Maximal value o the fidelity susceptibility Xmax as a function of system 
size L for two example three-body repulsions. On log-log scale the data points fit the 
linear predictions of finite scaling hypothesis. The slope coresponds to the ratio of two 
critical exponents 7/f. 



that the Hamiltonian (CQ) violates the universality hypothesis, or stated in another way 
that for the values of W considered, there is a marginal field with some overlap with 
the three body repulsion. 

What analyzed up to know is compatible with the presence of a line of conformally 
invariant fixed points connected by a marginal operator with some overlap with the 
three body operator. We can further check this scenario by considering the scaling 
of the entanglement entropy of a block. If the system is conformally invariant, the 
prefactor to the scaling provides an estimate of the central charges c of the underlying 
CFT. Indeed for systems with open boundary conditions the EE of a block of size / on 
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Table 1. Critical parameters of the studied system for different values of the three- 
body interactions W. Larger three-body repulsions stabilize dimer superfluid phase 
(larger critical tunneling J c ) and decreases critical exponent describing divergences of 
the fidelity susceptibility (7). Results are obtained with data collapsing method. 
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Figure 5. (color on-line) (upper panel) Scaling of the entanglement entropy (jlj 
as a function of the scaled block size for exemplary case W/\U\ = 30. QPT point 
corresponds to the vanishing of the oscillatory dependence of the parity of the block 
size, (bottom panel) Central charge of the system near the critical point for different 
values of three-body interactions. 



a finite lattice of length L takes the form §9\ EH SH 021 S31 IB] : 



S ± (x) 



6 



In 



sin 




+ s±(L) 




(4) 



including an alternating term s± (depending on the parity of the block sizes) arising due 
to the open boundary conditions as well as other corrections of order 0(l/L) including 
a boundary induced term. 

We calculate the EE scaling under open boundary conditions using the DMRG 
method for system size of L = 128 and consider its behaviour versus the first term of 
fllj). Deep in the single particle SF phase, there is a simple oscillating dependence on 
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the parity of the block size as seen for an exemplary value of VF/|C/| = 30, although 
the scaling dependence on In [sin(7rZ/L)] is approximately linear (Fig. [5]). On the other 
hand, the behaviour for J < J c is markedly different. The entropy does not apparently 
scale linearly with In [sin(7r//L)]. Significantly though, the oscillatory dependence on 
the parity of the block length is reversed compared to the case of the single particle 
SF. We therefore identify the change from PSF to SF (and thus the critical value 
J c ) with switching of the oscillating behaviour which also corresponds to the best 
possible linearization of the entropy dependence on In [sin(7r//L)]. This is true for all 
considered values of W. Importantly, the mentioned linearization also facilitates the 
direct estimation of the central charge for the QPT. The variation of the central charge 
with repulsion W is thus plotted in Fig. |5j For W — > oo we obtain a value very close 
to c = 3/2 as expected from previous results for the hard constraint model described 
by a U(l) x Z2 CFT. However, it is important to stress that having checked that the 
scaling of the EE has the same form than the one predicted in CFT is not sufficient 
to claim that the model is conformal invariant. There are several examples of systems 
known to exhibit similar scaling of the entanglement entropy to the one displayed by 
CFT even if they are not at all conformally invariant [4"5l 06] . In order to finally address 
the emergence of conformal invariance one should also check the scaling of the gap of 
the system and the form of the dispersion relation for small momenta so that we leave 
the issue about the presence of conformal symmetry at this QCP as an open question 
that deserves further studies of the model. Independentely on the presence of conformal 
symmetry, on moving away towards moderate values of the repulsion W, one observes 
(as in the critical exponent determined ealier) a growing dependence of the central 
charge on W again providing evidence towards non-universal behaviour in this model. 

4. Conclusions 

To conclude, in this article the one-dimensional attractive Bose-Hubbard model with 
three-body constrains was studied. The critical properties of the system were examined 
and the role of the strength of the on-site three-body repulsion was studied. The 
results presented here provide evidence that the tunable on-site three body repulsion in 
the attractive Bose-Hubbard model has both overlap with a relevant and a marginal 
operator, in the language of Renormalization Group theory, so that it both shift 
the critical point and influences the critical properties of the PSF-SF quantum phase 
transition. 
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